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FOREWORD 


rep°r^  rePresents  one  phase  of  a  general  study  of  the 

ARO  T^VVT  Wl«h  fluid  fl0WS  under  Contract  DA-31-124- 
^lkoff^ Reseawh  Office  -  Durham,  with  Dr.  Henry  R. 
veikoff  serving  as  principal  investigator. 
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ABSTRACT 


An  experimental  investigation  was  conducted  to  detect  the 
presence  of  electrically  induced  secondary  flows  inside  a  round 
pipe.  A  thin  concentric  wire  in  the  pipe  was  impressed  with  a 
high  positive  voltage  which  generated  a  corona  discharge,  pro¬ 
ducing  ions  and  a  radial  electric  field  in  the  pipe.  It  was 
postulated  that  secondary  flows  taking  the  form  of  vortex  flows 
inside  the  pipe  would  result. 

A  technique  of  utilizing  electrostatic  probes  inserted  into 
the  pipe  and  on  the  wall  provided  strong  evidence  that  secondary 
flows  were  generated  both  with  and  without  a  primary  axial  flow. 

At  zero  axial  flow  the  characteristic  length  of  the  secondary  flow 
was  on  the  order  of  two  pipe  diameters.  With  an  axial  flow,  when 
the  primary  and  secondary  flows  are  superimposed,  various  charac¬ 
teristic  lengths  were  observed. 
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STATEMENT  OF  THE  PROBLEM 

References  1-3  reported  significant  increases  in  friction  factor 
and  velocity  distribution  in  channel  flow  when  an  ion  current  was 
established  transverse  to  the  gas  flow.  As  an  ion  source  these  inves¬ 
tigations  utilized  a  corona  discharge  from  thin  wires  maintained  at  high 
voltage  with  respect  to  the  channel  walls;  ions  of  one  polarity  drifted 
out  of  the  corona  sheath  and  across  the  gas  stream  while  ions  of  the 
other  polarity  were  present  only  in  the  thin  corona  sheath  on  the  wires 
and  drifted  into  the  wire.  Guch  configurations  for  a  round  and  flat 
channel  are  shown  in  Fig.  1. 

Table  I  below  lists  some  of  the  results  of  the  effect  ot  the  ion 
current  on  the  main  channel  flow  friction  factor  for  laminar  flow  of  air. 
It  is  seen  that  first-order  effects  are  caused  by  the  ion-fluid  inter¬ 
action. 


Table  I 


Reference 

Channel 

Jo,  uA/inch2 

fi/f 

1 

1-|  inch  dia.  pipe 

0.94 

1.76 

2 

I4  inch  dia.  pipe 

1.77 

1.75 

3 

5/8"  x  5"  duct 

3.33 

1.8 

In  Table  I  above 

Jo  =  ion  current  density  at  the  channel  wall, 
f  =  friction  factor  for  no  current,  and 
fi  =  friction  factor  with  current. 

A  precise  theory  quantitatively  expressing  the  effect  of  the  trans¬ 
verse  ion  current  on  the  axial  channel  flow  has  not  yet  been  developed; 
in  fact,  the  mechanism  for  interaction  has  not  been  definitely  established. 
Reference  1  hypothesizes  that  the  corona  sheath  surrounding  the  wires 
directly  affects  the  mean  flow  through  viscous  action;  reference  2  postu¬ 
lated  the  existence  of  a  flow-induced  electric  field  which  results  in  an 
upstream-directed  body  force  on  the  mean  flow.  Reference  4  considered 
that  transverse  electric  forces  could  destabilize  the  mean  flow  and 
superimpose  electrically  induced  secondary  vortex  flows.  Such  secondary 
flows  would  then  provide  a  transport  mechanism  resembling  a  large-scale 
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Fig.  1  Typical  flow  channels  with  corona 
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turbulence  which  would  increase  the  axial  flew  friction  factor  and  change 
the  velocity  distribution. 


The  concern  of  the  present  investigation  is  to  experimentally  detect 
the  existence  of  electrically  induced  secondary  flows  by  means  of  electro¬ 
static  probes. 

An  electrostatic  probe  is  am  electrically  conducting  surface  which 
is  inserted  into  a  region  of  ionized  gas  and  through  collisions  of  ions 
on  the  probe  surface  draws  a  current  in  the  electric  circuit  connected 
to  the  probe.  If  different  values  of  potential  are  impressed  on  the 
probe  by  the  probe  circuit,  then  the  probe  current  will  be  affected  and 
consequently  a  probe  current-probe  voltage  curve  or  characteristic  is 
obtained.  If  the  probe  itself  does  not  significantly  change  the  properties 
of  the  ionized  gas,  then  from  the  probe  characteristic  the  number  density 
of  ions  and  the  potential  of  the  ionized  gas  at  the  location  of  the  probe 
can  be  calculated. 


DESCRIPTION  OF  THE  FLOW  CHANNEL  AND 
ELECTROSTATIC  PROBES 

The  flow  channel  used  in  the  present  report  was  1.5-inch  inside 
diameter  round  pipe  136  inches  long.  A  0.004-inch  diameter  steel  wire 
w as  stretched  taut  along  the  centerline  and  insulated  from  the  pipe  walls. 
At  a  point  115  inches  downstream  from  one  end  of  the  pipe  a  l/8-inch 
diameter  brass  rod  probe  was  inserted  1/8  inch  into  the  pipe  along  a 
radius.  The  sides  of  this  rod  were  covered  with  a  covering  of  l/l6-inch 
thick  nylon,  so  only  the  end  of  the  rod  surface  contacted  the  gas  in  the 
pipe.  One-half  inch  further  downstream  was  a  similarly  placed  rod  probe 
identical  to  the  one  above.  The  end  surfaces  of  these  two  rods  served 
as  electrostatic  probe  surfaces  and  will  hereafter  be  referred  to  as  the 
rod  probes.  A  sketch  of  the  rod  probes  is  shown  in  Fig.  2, 

A  section  of  the  pipe  beginning  approximately  120  inches  downstream 
from  the  end  of  the  pipe  was  constructed  of  a  series  of  thirty  3/32-inch 
wide  rings,  each  separated  electrically  by  an  0.004-inch  thick  intervening 
sheet  of  vinyl  plastic.  Each  ring  had  a  pin  imbedded  in  its  outside 
diameter  to  serve  as  an  electrical  connector,  and  the  entire  set  of  rings 
were  cemented  together  and  positioned  in  line  with  the  rest  of  the  pipe 
(see  Fig.  2).  There  were  30  rings  in  all,  numbered  consecutively  from 
the  downstream  ring  to  the  upstream  ring.  The  inside  surface  of  any 
ring  of  course  was  exposed  to  the  gas  in  the  pipe  and,  therefore,  could 
be  used  as  an  electrostatic  probe  whose  position  is  at  the  radius  R  of 
the  pipe. 

Air  was  supplied  to  the  flow  channel  from  the  laboratory  compressed 
air  supply,  as  shown  in  Fig.  3»  The  downstream  end  of  the  pipe  exhausted 
into  the  atmosphere.  The  central  corona  wire  was  connected  to  a  high- 
voltage  power  supply;  the  pipe  itself  (aluminum)  was  electrically 
grounded. 
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ELECTROSTATIC  EQUATIONS 

The  electric  field  equations  are  assumed  to  one -dimensional;  i.e., 
no  component  of  field  strength  exists  in  the  axial  direction.  Using 
Maxwell's  equation 


— >  n 
E  =  | 


(1) 


or 


1  d(Er)  .  n  . 
r  dr  e 


(2) 


The  local  current  density  J  is  given  by 


J  =  nkE  , 


(3) 


where 


n  is  the  number  of  density  of  positive  ions, 
k  is  the  ion  mobility,  and 
E  is  the  local  value  of  field  strength. 

From  current  continuity 

Jr  =  JoR  (4) 

where  subscript  o  represents  values  at  the  pipe  wall,  r  *  R.  Equa¬ 
tion  (2)  then  reduces  to 

(Er)d(Er)  =  ^  r  •  (5) 

Integrating, 

(Er)2  =  r2  +  constant.  (6) 
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A  space-charged  limited  condition  is  assumed  at  the  edge  of  the  corona 
sheath  at  r  =  ri.  Therefore,  E  =  0  at  r  =  n,  Equation  (6)  yields 


(Er)2  =  Cr2  -  n2]  ;  (7) 

but  the  corona  sheath  radius  n  is  small  compared  to  the  radius  r  for 
most  of  the  pipe,  so  Eq.  (7)  is  written  approximately  as 


Hie  ion  density  n,  using  Eq,  (3),  (4),  and  (8),  can  be  written 


Equation^  9)  can  be  averaged  over  the  pipe  volume  to  find  an  average  ion 
density  n: 


n  » 


(10) 


Following  reference  2  a  nondimens ional  group,  called  a  charge 
number,  is  defined  as 


N 


T*Jocl»  R* 
.  Rk  J  Mk 


For  the  flow  channel  in  this  report  the  charge  number  as  a  function 
of  current  density  Jo  on  the  pipe  wall  is  given  by  Table  II,  where  a 
positive  ion  mobility  of  0,216  inch2  volt"1  sec"1  was  assumed,5 


Table  II 


Charge  Number 

N 

Current  Density 

H  A/inch2 

2 

O.OO963 

5 

0.0602 

10 

0.241 

15 

0.542 
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The  charge  number  N  was  found  in  Ref.  2  and  3  to  be  a  convenient 
parameter  for  correlating  the  ion  current  effect  on  the  friction  factor. 
For  example ,  in  the  present  case  the  static  pressure  drop  between  two 
pressure  taps  at  45  inches  and  127  inches  from  the  upstream  end  of  the 
channel,  and  the  friction  factor  were  calculated.  In  laminar  flow,  the 
friction  factor  at  a  constant  charge  number  is  seen  to  possess  a  laminar 
characteristic,  as  shown  in  Fig.  4.  Also,  the  group  fi  Rey  is  found  to 
be  linear  with  the  charge  number,  as  shown  in  Fig.  5.  These  results 
agree  with  the  findings  reported  in  Ref,  2  and  3. 


MEASUREMENT  OF  ION  DENSITY 

The  probe  circuit  shown  in  Fig.  6  was  first  connected  to  the  up¬ 
stream  rod  probe.  The  potential  of  the  probe  was  varied  by  adjusting 
the  variable  resistors  shown  in  Fig.  6.  For  no  corona  discharge  in  the 
pipe,  the  measured  current  due  to  leakage  was  very  small  and  on  the  order 
of  0.0034  A  at  a  probe  potential  of  400  V  relative  to  the  pipe  wall.  A 
corona  discharge  was  then  generated  on  the  axial  corona  wire  and  the 
corona  current  set  a  value  corresponding  to  a  current  density  of  0.43  4A 
per  square  inch  on  the  pipe  wall.  From  Eq.  (ll)  this  corresponds  to  a 
charge  number  of  13.05,  well  within  the  region  of  significant  ion-fluid 
interaction  as  seen  from  Fig.  5.  With  this  value  of  charge  number  the 
probe  characteristic  was  obtained  for  an  axial  flow  Reynolds  number  of 
0,  1270,  and  3400  and  is  shown  in  Fig.  7.  During  this  test  the  charge 
number  varied  less  than  1 %  because  of  slight  drift  in  the  value  of  the 
corona  current. 

At  high  positive  values  of  probe  voltage  the  probe  repels  the 
positive  charges  and  the  probe  current  is  low.  Under  this  condition 
the  gas  near  the  probe  surface  is  partially  depleted  of  positive  ions 
as  compared  to  the  undisturbed  ion  density.  At  high  negative  values  of 
probe  voltage,  the  probe  begins  to  attract  all  ions  in  the  vicinity  of 
the  probe  and  the  current  tends  to  level  off  to  a  high  value. 

If  it  is  assumed  that  all  ions  which  strike  the  probe  are  measured 
as  current,  that  the  ions  possess  a  translational  temperature  very  near 
the  gas  temperature,  and  that  the  velocity  distribution  is  Maxwell- 
Boltzmann,  then  the  ion  density  can  be  calculated  from  Fig.  7»5  At  a 
point  on  the  probe  characteristic  corresponding  to  the  knee  of  the  curve, 
the  probe  is  at  a  potential  equal  to  the  space  potential  of  the  ion  gas 
at  that  point.  Under  this  condition  the  probe  neither  attracts  nor 
repels  ions  and  the  probe  current  Jp  is  related  to  the  local  ion  density 
by  the  equation  given  in  Ref.  5>  page  140. 


n  =  4jp 


(12) 
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Pig.  5  Friction  Factor  Times  Reynolds  Number  Versus  Charge  Number 


Electrostatic  probe  circuit 


Probe  Current 


where 


m  is  the  mass  of  the  ion, 

T  is  temperature,  and 
K  is  Boltzmann's  constant. 

Since  the  knee  of  the  curve  is  not  well  defined,  the  probe  current 
density  can  only  be  determined  within  about  a  factor  of  10.  For  Fig.  7 
the  knee  is  at  approximately  a  probe  voltage  of  250. 

The  three  sets  of  points  on  Fig.  7  for  different  axial  flow  Reynolds 
numbers  apparently  fall  on  the  same  curve.  There  is,  therefore,  no 
indication  that  the  mean  axial  flow  of  air  has  any  effect  on  the  local 
field  strength  or  ion  density  at  the  point  of  the  probe.  This  conclusion 
is  further  supported  by  Fig.  8  which  shows  the  change  in  probe  current 
with  Reynolds  number  for  a  fixed  probe  voltage  of  240.  No  significant 
change  in  probe  current  is  observed  as  the  Reynolds  number  is  varied, 
except  for  a  very  small  increase  near  a  Reynolds  number  of  zero. 

Ring  number  15  was  also  used  as  an  electrostatic  probe,  with  all 
'  the  other  rings  grounded.  In  this  case,  as  the  probe  voltage  was 

changed  relative  to  ground,  the  leakage  current  was  measured  as  a  func¬ 
tion  of  probe  voltage  (see  Fig.  9)»  A  corona  current  was  generated  in 
the  pipe  corresponding  to  a  charge  number  of  18.3  and  the  probe  current 
measured  as  a  function  of  the  probe  voltage.  The  measured  probe  currents 
were  then  corrected  for  leakage  (see  Fig.  9),  and  the  net  probe  current 
plotted  against  the  probe  voltage  as  shown  in  Fig.  10.  In  this  case  the 
probe  is  at  the  outer  radius  R  and  the  space  potential  should  therefore 
be  a  a  value  of  zero.  The  zero  probe  voltage  point  on  Fig.  10  is  with 
the  portion  of  the  curve  constituting  the  knee,  as  expected.  An  almost 
identical  curve  was  obtained  for  three  Reynolds  numbers  of  0,  1260,  and 
3400,  again  indicating,  as  with  the  rod  probe,  that  the  ion  density  is 
not  measurably  changed  by  the  mean  air  flow.  This  is  similarly  shown  by 
a  plot  of  the  probe  current  versus  Reynolds  number  at  a  constant  probe 
voltage  (Fig.  11).  No  large  changes  in  Jp  are  observed  as  the  Reynolds 
number  is  varied  up  to  a  value  of  5000. 

The  ion  density  was  calculated  from  the  probe  characteristics  and 
Eq.  (12)  and  was  compared  with  the  theoretical  value  from  Eq.  (9).  The 
results  for  the  rod  and  ring  probes  are  shown  in  Table  III.  A  tempera¬ 
ture  of  70°F  and  a  mobility  of  0.216  inch2  volt-1  sec-1  were  used  in  the 
calculations.  This  value  is  representative  of  the  mobility  of  the  pos¬ 
itive  ion  in  air  as  given  by  Ref.  5,  page  38. 

The  measured  ion  density  values  in  Table  II  indicate  reasonable 
agreement  with  the  values  predicted  by  Eq.  (9).  In  particular  it  is 
noted  that  the  level  of  ion  density  is  very  low  compared  to  the  number 
density  of  neutral  molecules;  i.e.,  on  the  order  of  one  to  10. 11 
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Fig.  10  Probe  Characteristic  for  Ring  Probe,  N 
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Probe  Current  Versus  Reynolds  Number,  ring  probe 


Table  III.  Ion  Density 


Ion  Density,  Ions/inch3 

N 

Current  Density  J0 
(iA/inch2 

Eq.  (9) 

Rod  Probe 

Ring  Probe 

18.3 

.80 

6.6  x  109 

— 

1  x  109 

13.05 

.41 

5.6  x  109 

2  x  109 

— 

MEASUREMENT  OF  MOVING  DISTURBANCES 
1.  Rod  Probe 

The  radial  current  density  J  results  in  a  body  force  per  unit 
volume  in  the  radial  direction.  The  equations  of  motion  admit  a  solu¬ 
tion  such  that  no  radial  velocity  component  exists  and  a  steady  radial 
pressure  gradient  J/k  is  present.  It  is  hypothesized,  however,  that 
the  primary  flow  is  unstable,  and  a  secondary  flow  consisting  of  a 
toroidal  vortex  pattern  (as  in  Fig.  12)  will  be  superimposed  on  the 
primary  flow.  The  secondary  and  the  primary  flows  will  in  general 
combine  in  a  complicated  way. 

A  method  for  detecting  the  presence  of  the  secondary  flow 
results  from  the  following  reasoning: 

When  the  primary  flow  is  zero  axial  motion,  the  only  motion 
present  is  the  secondary  flow.  This  type  of  flow  has  an  analogy  with 
the  case  of  fluid  between  concentric  cylinders  with  the  inner  cylinder 
rotating.  The  secondary  flow  in  that  case6  was  found  to  be  toroidal 
vortices  with  a  diameter  equal  to  twice  the  width  of  the  annulus.  If 
in  the  present  case  a  similar  pattern  of  vortices  is  generated  at  zero 
primary  flow,  then  an  observer  moving  along  the  pipe  at  velocity  v 
would  encounter  variable  velocities  with  a  periodicity  Vv>  where  A  is 
the  wavelength.  It  is  now  assumed  that  if  the  primary  flow  is  not  zero, 
the  combination  of  the  primary  and  secondary  flow  will  be  such  that  a 
periodicity  in  the  resultant  flow  will  be  maintained.  Under  this  condi¬ 
tion  a  stationary  observer  in  the  pipe  will  experience  a  periodic  dis¬ 
turbance  in  velocity  with  a  period  Vv>  where  v  is  now  the  velocity  of 
the  resultant  flow  and  }\  is  the  wavelength.  This  is  akin  to  saying  that 
the  primary  flow  carries  the  secondary  flow  downstream. 

The  velocities  associated  with  the  secondary  flow  are  expected 
to  be  very  small  and  not  detectable  by  ordinary  means.  However,  by 
looking  at  the  secondary  flow  pattern  in  Fig.  12,  it  is  seen  that  as 
the  pattern  is  moved  past  a  stationary  observer,  the  velocity  will 
alternately  add  and  detract  from  the  radial  drift  velocity  of  the  ions. 
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Fig.  12  Sketch  of  Induced  Vortex  Flow 


Such  an  occurrence  would  impart  a  periodicity  to  the  local  ion  density 
and  current  at  the  observer's  position. 

The  two  rod  probes  (Fig.  2)  were  therefore  connected  to  the 
terminals  of  an  oscilloscope.  If  no  secondary  flow  were  present  it 
would  be  expected  that  no  pattern  would  be  detected  on  the  scope; 
however,  if  secondary  flows  were  present,  from  the  previous  reasoning  a 
periodic  disturbance  would  bo  detected.  Figures  13-18  show  the  resulting 
oscilloscope  traces  for  various  Reynolds  numbers  from  0  to  1780  for  a 
charge  number  12.75*  Photographs  which  show  two  traces  merely  indicate 
that  two  sweeps  were  made  on  the  same  film.  From  these  traces  the 
following  observations  are  made:  For  zero  flow,  a  slight  periodic  dis¬ 
turbance  is  measured;  as  the  Reynolds  number  is  increased,  the  amplitude 
of  the  distrubances  increases  to  maximum  values  in  the  range  of  600  to 
900  Reynolds  number  and  then  decrease  again  to  small  relatively  constant 
amplitudes  above  a  Reynolds  number  of  1330. 

The  time-dependency  of  the  traces  reveals  that  the  disturb¬ 
ances  are  not  regular;  however,  certain  recurring  wavelengths  are 
observed.  It  was  anticipated  that  the  moving  disturbances  would  have  a 
wavelength  equal  to  the  pipe  diameter  and  move  with  the  bulk  fluid  at  a 
velocity  equal  to  the  mean  velocity  of  the  fluid.  Therefore,  the  fol¬ 
lowing  length  was  calculated 


(Pipe  Diameter)  (Sweep  Speed) 
(Mean  Velocity) 


and  appears  as  a  length  indicated  on  each  photograph.  This  distance  is 
the  horizontal  length  on  the  scope  screen  which  corresponds  to  moving 
disturbances  with  a  characteristic  length  equal  to  one  pipe  diameter. 

From  the  scope  traces  it  is  seen  that  the  distance  between  sucessive 
relative  maxima  on  the  traces  frequently  does  approximately  equal  the 
pipe  diameter.  This  fact  indicates  that  disturbances  with  a  wavelength 
of  one  pipe  diameter  are  moving  down  the  pipe.  It  is  also  noticed  that 
waves  with  longer  and  shorter  lengths  are  also  present. 

Above  a  Reynolds  number  of  about  1150  the  disturbance  pattern 
is  seen  to  begin  to  posses  a  high-frequency  con5>onent,  which  also  is 
accompanied  by  a  decrease  and  finally  a  leveling  off  of  the  amplitude. 

It  may  be  that  in  this  range  of  flow  the  secondary  flow  and  the  primary 
flow  combine  in  a  different  way,  and  the  relatively  more  regular  pattern 
at  low  Reynolds  numbers  is  lost.  It  may  also  be  suspected  that  transition 
to  true  turbulence  begins  in  this  range;  however,  Fig.  4  indicates  that 
the  friction  factor  exhibits  a  laminar  characteristic  up  to  a  Reynolds 
number  of  more  than  2000.  It  is  also  noted  that  the  friction  factor 
curve  of  Fig.  4  did  not  seem  to  experience  any  change  in  trends  at  any 
Reynolds  number  below  the  transition  Reynolds  number  at  about  2500. 
Another  possibility  is  that  beginning  at  Reynolds  numbers  of  about  1150, 
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Fir.  13  Fluctuations  Between  Rod  Probes 
=  12.75;  Sweep  =  0.1  sec/cm;  Amplitude  *  1  mv/cm 
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Fir.  14  Fluctuations  Between  Hod  Probes 
N  =  12.75;  Sweep  =  0.1  sec/cm;  Amplitude  *  2  mv/cm 


Fig.  15  Fluctuations  Between  Rod  Probes 
12.75;  Sweep  ■  0.1  sec/cm;  Amplitude  =  2  mv/cm 


Amplitude  »  2  mv/cm 


Fig.  16  Fluctuations  Between  Rod  Probes 
N  ■  12.75;  Sweep  «  0.1  sec/cm 
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Fig.  18  Fluctuations  Between  Rod  Probes 


the  presence  of  the  probes  induces  a  true  turbulence  in  the  region  near 
the  probes  which  is  reflected  in  a  high-frequency  component  to  the 
observed  disturbances;  yet  transition  to  turbulence  in  the  major  portion 
of  the  pipe  has  not  occurred. 

The  last  photograph  on  Fig.  18  shows  the  trace  for  no  flow 
and  no  corona  current  to  give,  as  a  matter  of  reference,  the  background 
noise  level. 

Since  the  disturbance  patterns  were  observed  to  be  combinations 
of  many  wavelengths,  a  harmonic  analysis  was  performed  with  the  aid  of 
a  computer.  Figures  19  and  20  show  the  traces  which  were  analyzed  and 
the  results  of  che  analysis.  The  wavelengths  listed  correspond  to  the 
length  of  a  disturbance  moving  down  the  pipe  at  a  velocity  equal  to  the 
mean  velocity  of  the  axial  flow.  The  relative  magnitudes  of  certain 
wavelengths  present  in  the  trace  are  listed.  The  results  indicate  that 
the  main  disturbances  have  a  characteristic  length  of  about  2,  3>  or  4 
times  the  pipe  diameter. 

Of  course  since  the  disturbances  are  measured  by  two  stationary 
probes  a  fixed  distance  of  |--inch  part,  the  relative  magnitudes  of 
various  wavelengths  will  be  distorted  when  measured  by  the  probes. 

2.  Detection  of  a  Critical  Charge  Number 

If  the  occurrence  of  secondary  flow  is  truly  caused  by  a 
destabilization  of  the  primary  flow,  a  critical  value  of  the  destabi¬ 
lizing  electric  parameter  would  be  expected  to  exist,  below  which  the 
primary  flow  is  stable.  The  ion  current  effect  on  the  friction  factor, 
and  the  disturbances  measured  on  the  probes  should  therefore  theoretically 
not  appear  until  a  finite  value  of  current  density  is  exceeded.  The 
critical  current  density  may  also  depend  on  the  Reynolds  number. 

An  attempt  to  detect  a  critical  charge  number  was  performed  as 

follows ; 

A  constant  axial  flow  Reynolds  number  of  1150  was  obtained  in 
the  pipe.  With  the  two  rod  probes  connected  to  the  oscilloscope,  the 
disturbance  pattern  was  observed  as  the  ion  current  was  increased.  At 
the  same  time,  the  static  pressure  drop  in  the  pipe  between  points  1*5 
inches  and  127  inches  from  the  upstream  end  were  measured.  The  resulting 
oscilloscope  traces  are  shown  in  Fig.  21.  A  very  slight  p*'v  dicity  is 
seen  to  appear  at  the  fourth  trace,  at  a  charge  number  of  3.52;  at  higher 
charge  numbers  the  distrubance  pattern  is  more  obvious.  Below  N  =  3»52 
a  wave  pattern  is  not  observed,  however,  it  is  realized  that  the  waves 
may  be  present  at  small  amplitude.  The  existence  of  the  relatively 
large  amplitude  of  high-frequency  corona  oscillations  in  this  range  makes 
the  detection  of  onset  difficult. 

Figure  22  chows  the  friction  factors  measured  as  a  function 
of  charge  number  in  this  range.  A  definite  increase  in  the  friction 
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Fig.  19  Harmonic  Analysis  of  Oscilloscope  Trace 
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factor  is  observed,  and  the  trend  is  to  be  a  linear  fr.u  cion  of  the 
charge  number.  If  the  line  at  the  higher  charge  numbers  is  extra¬ 
polated  to  the  value  of  the  friction  factor  at  zero  charge  number,  the 
charge  number  at  this  point  represents  the  critical  charge  number,  or 
onset  of  electric  effect  on  the  friction  factor.  The  dashed  line  on 
Fig.  22  shows  this  point  to  be  at  a  charge  number  of  1.8,  corresponding 
to  a  current  density  of  1.52  uA/inch2  of  pipe. 

The  detection  of  a  critical  current  is  a  difficult  experimental 
task.  At  very  low  current  densities  in  the  pipe  the  current  distribution 
along  the  pipe  tends  to  become  nonuniform  because  of  slight  changes  in 
geometry  along  the  pipe.  The  measurement  of  the  friction  factor  ever 
a  relatively  long  length  of  pipe  may  then  conceivably  include  sections 
of  pipe  which  are  above  critical  current  and  sections  where  the  current 
is  below  critical.  The  above  method  of  using  probes  to  detect  secondary- 
flow  is  similar  in  principle  to  a  method  used  by  Donnelly7  to  detect 
vortices  between  concentric  cylinders.  In  the  present  case,  however, 
the  probe  method  is  made  difficult  because  the  ion  drift  velocity  is  very 
large  compared  to  the  velocities  characteristic  of  secondary  flows. 

3.  Ring  Probes 

Disturbances  were  also  measured  by  having  two  of  the  ring 
probes  connected  to  the  oscilloscope.  The  extreme  downstream  ring 
(ring  no.  l)  and  the  eight  ring  were  used  as  probes;  the  probe  separa¬ 
tion  distance  was  therefore  0.68  inch.  Figures  23  and  24  show  the 
observed  disturbance  patterns  at  a  constant  charge  number  of  18.3  for 
several  different  Reynolds  numbers.  Again  at  zero  Reynolds  number  a 
low-frequency  oscillation  is  seen  to  be  present.  At  finite  Reynolds 
numbers  the  disturbance  patterns  again  reveal  characteristic  wavelengths. 
The  length  indicated  on  the  photographs  indicates  the  distance  corre¬ 
sponding  to  a  wavelength  of  one  pipe  diameter,  based  on  the  mean  axial 
flow  velocity  as  before.  Several  traces,  such  as  the  next  to  the  last 
one  on  Fig.  23  and  24  indicate  that  the  main  disturbance  has  a  wave¬ 
length  approximately  twice  the  pipe  diameter.  In  any  case,  the  results 
do  indicate  that  characteristic  wavelengths  are  indeed  present. 

The  amplitudes  of  the  disturbances  were  a  maximum  for  the 
traces  at  Rey  =  900,  and  decreased  at  higher  Reynolds  numbers,  in 
agreement  with  the  trend  observed  with  the  rod  probes.  The  amount  of 
decrease  with  the  ring  probes  was  less,  however.  Another  difference 
in  the  observed  traces  with  the  two  sets  of  probes  is  the  type  of  dis¬ 
turbances  at  high  Reynolds  numbers.  Beginning  at  around  Rey  =  1300, 
the  disturbances  with  the  rod  probes  began  to  contain  a  high-frequency 
component;  this  phenomenon  was  not  observed  with  the  ring  probes.  It 
is  not  known  at  this  time  whether  this  difference  is  caused  by  the 
disturbing  effect  of  the  rod  probes  themselves,  or  the  true  disturbance 
patterns  at  different  pipe  radii. 

Figure  25  shows  similar  traces,  but  at  a  constant  Reynolds 
number  and  increasing  charge  number.  The  lowest  charge  number  indicated 
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Fir.  23  Fluctuations  Eetween  Ring  Probes 
K  =  18.3;  Amplitude  *  5  mv/cm 
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Fig.  24  Fluctuations  Between  Ring  Probes 
N  =  18.3;  Amplitude  ■  5  mv/cm 
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a  characteristic  wavelength  of  about  four  diameters;  at  higher  charge 
numbers  wavelengths  of  the  order  of  one  diameter  are  discemable  in  the 
traces. 


MEASUREMENT  OF  STATIONARY  DISTURBANCES 

Toroidal-shaped,  electrically  induced  vortices  are  postulated  to 
be  generated  in  the  pipe  at  a  Reynolds  number  of  zero.  When  the  primary 
flow  consists  of  an  axial  motion  this  secondary  flow  would  combine  with 
the  Poiseuille  primary  flow  in  a  complicated  way.  It  would  be  advan¬ 
tageous  then  to  be  able  to  measure  the  secondary  flow  when  the  Reynolds 
number  is  zero.  To  achieve  this  end,  part  of  the  current  from  a  ring 
was  passed  through  a  resistor  before  being  grounded,  and  the  resulting 
voltage  drop  across  the  resistor  was  placed  across  the  terminals  of  an 
oscilloscope,  A  mechanically  driven  rotary  switch  successively  connected 
each  ring  to  the  resistor,  so  that  the  net  effect  over  one  revolution 
of  the  switch  was  to  indicate  on  the  oscilloscope  the  relative  variation 
of  current  in  the  rings  (see  rotary  switch  circuit,  Fig.  26).  The 
rotary  switch  turned  one  revolution  every  1,35  seconds,  A  variable 
amount  of  current,  due  to  secondary  flow,  would  hopefully  be  reflected 
in  a  disturbance  pattern  on  the  oscilloscope. 

Figures  27-28  show  the  resulting  oscilloscope  traces.  Within  the 
range  of  the  moving  probe,  a  peak  on  the  traces  is  observed,  suggestive 
of  a  characteristic  wavelength  on  the  order  of  two  pipe  diameters,  TMs 
measurement  has  the  drawback  that  the  scanned  portion  of  the  pipe  is 
less  than  one  wavelength.  If  an  axial  flow  is  present  when  the  rotating 
switch  is  operating,  the  resulting  traces  are  as  shown  in  Fig.  29,  This 
case  is  a  combination  of  a  moving  and  stationary  probe,  where  the  charac¬ 
teristic  velocity  would  be  the  sum  of  the  mean  air  velocity  and  the 
effective  velocity  of  the  moving  ring  probe. 

Since  the  rotary  switch  was  limited  to  contacting  23  variable  ring 
probes,  a  longer  portion  of  the  pipe  could  be  scanned  by  the  following 
method.  Every  odd-numbered  ring  was  connected  to  the  rotary  switch 
while  every  even-numbered  ring  was  grounded.  In  this  way  a  total  length 
of  2,93  inches  of  the  pipe  was  scanned  by  the  moving  probe.  Figure  30 
shows  the  traces  observed  for  several  charg  numbers.  When  the  rotating 
contact  is  between  two  contacts,  the  oscilloscope  beam  moves  to  the  zero 
position  (upward  in  the  photographs);  the  traces,  therefore,  appear  as 
a  nearly  vertical  lines,  the  terminals  of  which  represent  the  measured 
quantity.  The  traces  on  Fig.  30  show  that  almost  one  complete  wave¬ 
length  of  about  2.8  inches,  or  very  nearly  two  pipe  diameters,  is  present. 
Figure  31  shows  a  similar  series  of  photographs  but  in  this  case  the 
oscilloscope  was  fixed  with  a  low-persistance  screen;  as  a  result  the 
only  visible  trace  appears  when  the  electrcr.  beam  dwells  for  an  instant 
before  moving  to  the  zero  position.  These  points  show  some  scatter, 
but  again  almost  a  complete  wavelength  corresponding  to  about  two  pipe 
diameters  is  observed. 
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Fig.  26  Rotary  Switch  Circuit 


Fig.  27  Fluctuations  on  Moving  Ring  Probe 
Rey  ■  0;  Sweep  ■  .5  sec/cm 
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>8  Fluctuations  on  Moving  Ring  Probe 
Rey  =  0;  Sweep  B  .5  sec/cm 


Fig.  29  Fluctuations  on  Moving  Ring  Probe 
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SUMMARY 


The  electrostatic  probe  technique  used  in  this  investigation  showed 
that  local  disturbances  in  the  electric  field  were  present  in  the  pipe. 
These  disturbances  were  postulated  to  be  caused  by  coupling  of  an 
electrically  induced  vortex  motion  with  the  ion  current.  In  axial  flow, 
such  induced  vortex  motion  would  be  expected  to  considerably  affect  the 
momentum  transport  within  the  pipe,  in  agreement  with  the  experimental 
observations  of  the  increase  in  friction  factor.  From  the  wavelength 
of  the  disturbances,  there  were  indications  that  the  secondary  flows 
had  a  characteristic  length  of  about  two  pipe  diameters  when  no  axial 
flow  was  present,  and  that  several  characteristic  lengths  were  super¬ 
imposed  when  there  was  axial  flow. 

A  determination  of  the  local  ion  density  at  a  point  within  and  at 
the  wall  of  the  pipe  using  the  probes  gave  substantial  agreement  with 
the  theoretical  one -dimensional  equations.  No  effect  on  the  local  ion 
density  was  measurable  as  the  axial  flow  Reynolds  number  was  increased 
from  0  up  to  a  value  of  over  5000. 

The  results  of  the  present  investigation  give  support  to  the  concept 
of  an  electric  destabilization  of  flow  in  a  pipe  as  proposed  in  Ref.  4. 
Since  the  interaction  of  the  secondary  flow  with  the  primary  flow  is 
quite  large,  it  is  recommended  that  the  analysis  be  extended  to  cover 
additional  geometrical  and  electrical  configurations. 
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